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Phase segregation for binary mixtures of Bose-Einstein 

Condensates 

M. Goldman* B. Merletl 


Abstract 

We study the strong segregation limit for mixtures of Bose-Einstein condensates modeled 
by a Gross-Pitaievskii functional. Our first main result is that in presence of a trapping 
potential, for different intracomponent strengths, the Thomas-Fermi limit is sufficient to 
determine the shape of the minimizers. Our second main result is that for asymptotically 
equal intracomponent strengths, one needs to go to the next order. The relevant limit is 
a weighted isoperimetric problem. We then study the minimizers of this limit problem, 
proving radial symmetry or symmetry breaking for different values of the parameters. We 
finally show that in the absence of a confining potential, even for non-equal intracomponent 
strengths, one needs to study a related isoperimetric problem to gain information about the 
shape of the minimizers. 


1 Introduction 

In this paper, we investigate the asymptotic behavior as e goes to zero of the Gross-Pitaievskii 
functional 


Feiv) = £ [ |Vt/iP- b |Vr/2|^ + 1 / lvt + ^rit + Kr]lrj‘l + {T]l + r]j)V, (1) 

Jr2 e Jr2 2 2 

under the mass constraint 



This functional arises in the study of two component Bose-Einstein condensates. It has been 
widely studied, both in the physical and mathematical literature (see [25] and the references 
therein or the book |2|). The potential E is a trapping potential. For simplicity, we only consider 
here the harmonic potential V = \x'\^. The constant g, measures the asymmetry between the 
intracomponent repulsive strengths of each component and K represents the intercomponent 
repulsive strength. Without loss of generality, we will take here g > 1. The case K < y/g, where 
mixing of the two condensates occurs has recently been well understood in |3]. On the other 
hand, the case K > y/g, where it is expected both experimentally [331 [SS], numerically [331 [37] 
and theoretically [sgiiniiiiiiiiiiis] that segregation occurs, is maybe not yet so well understood. 
In the symmetric case g = 1, it has been proved in [11E5] that, as expected from the physics 
literature, the Thomas-Fermi limit i.e. the limit of sF^ only imposes segregation but does not 
give any information about the actual shape of the minimizers. To gain such information, one 
has to go to the next order since — min^ converges in the sense of T—convergence |14j 
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to a weighted isoperimetric problem. In this paper, we focus on the asymmetric case g > 1 and 
show that the situation is radically different. 

Let us introduce the Thomas-Fermi energy: 



1 

2 


Pi + ^P2 + Kpip2 + {pi + P2)V 


Our first main theorem is the following: 


( 3 ) 


Theorem 1.1. For every ai,a 2 ,g,K > 0, with K > ^^g > 1, there exists a unique minimizer 
P° = (Pi,P 2 ) of ([3]) under the volume constraints 



and 



( 4 ) 


This minimizer p^ is radially symmetric (see Lemma \2.1\ for explicit formulas for p^ and Eq = 
E{p^)). 

Moreover, we have the following stability result: there exists C > 0 (which depends only on 
01,02 and g) such that if p satisfies the constraints (H]) thei^ 


\\p-p^\\i<CVE{p)-Eo. 


( 5 ) 


Let us point out that the radial symmetry of the minimizer of the Thomas-Fermi energy was 
not completely expected (see [HI [32] or the discussion in [H Sec. 1.3.4]). As a consequence of 
this stability result, we prove that minimizers of E^ converge to V^)- 

Theorem 1.2. There exists C > 0 such that for e G (0,1], any minimizer of E^ under the 
constraints m satisfies 

< ( 6 ) 

2 

This theorem establishes that in the non-symmetric case, the Thomas-Fermi limit already 
provides full information on the limiting behavior of the minimizers of ([T]). It is quite surprising 
that even without using isoperimetric effects, we are able to obtain strong convergence of the 
minimizers in the form of Q. Let us point out that the idea of using stability inequalities such 
as ([5]) to get (quantitative) convergence results is far from new (see for instance [13 [12]). One 
crucial point which explains the difference between the asymmetric case and the symmetric one 
is that here, there is a gap between the two Thomas-Fermi profiles p^ and p® the sense that 
there exists ro > 0 such that supppj C i?(0, ro) and supppg C \ B{0,ro) with 

inf p?(r) > supp 2 (’")- 

r<rQ r>ro 



We then study the crossover case where g = 1 + for some .^ > 0. Let be the minimizer 
of 

Ge{r])=£[ + 

J'R? £ ^ 

under the mass constraint 



^we denote by || • ||p the U’ norm 
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where a = ai + 02 - It is well known [261 M that 'q'l converges when e goes to zero to the 
Thomas-Fermi profile 

-p={R^-V)+ 

where R is chosen such that f ^2 p = a. Building on results obtained in the case g = 1 |25l 0] , 
we get 


Theorem 1.3. For K > the functional {F^ — F— converges as e goes to zero for the 

strong topology to 


G^iui,U2) = 


I + ^ [ p^ 

IdE Je<^ 

+00 


ifui=pxE, U2 =pxe<^ and / p = ai, 

JE 

otherwise, 


where ax > 0 is defined by the one dimensional optimal transition problem 

= inf + |r/ 2 |^ + ^ - 1 )^ + (iF - l)r7i?72 : • 


We study the minimizers of the limiting functional and prove our second main result. 
Theorem 1.4. The following holds: 

• there exists oq S (0, a/2] such that for every ai € (ao, a — oq) there exists such that 

the minimizer of is not radially symmetric for ^ , 

• for every ai € (0,a), there exists such that the minimizer of is the centered ball for 

The regime g = 1 + corresponds to the numerical simulation of [32] , In that paper, the 
observed numerical results (droplets) fit with the first point of Theorem 11.41 (see in particular 
|32l Fig l.c, Fig. 3.a]). The first part of Theorem 11.41 is a consequence of a symmetry breaking 
result from |3]. The second part follows from a combination of two results. The first is a stability 
result for the functional fP ': 

Proposition 1.5. For every a G (0,a), there exists C = C{a) > 0 such that for every measur¬ 
able set E C with f^p = a, we have, 



where r is such that Jg p = a. 

The second is an estimate on the potential instability of the ball for the weighted isoperimetric 
problem. 

Proposition 1.6. For every a G (0, a), there exists c = c(a) > 0 such that for every set E C 
with locally finite perimeter and with f^p = a, there holds 

f fG. f fG > -c([ p)\ (7) 

JdE JdBr KJEABr / 

where r is such that p = a. 


3 


The rigidity result given by Theorem ll.4l is similar in spirit to several rigidity results obtained 
for variants of isoperimetric problems (see |29[I21[[2^ for instance). However, the peculiar aspect 
here is that this rigidity does not come from the isoperimetric term but rather from the volume 
term. Nevertheless, the proof of ([7]) bounding the instability of the ball follows the strategy of 
m (see also mm) to prove quantitative stability estimates for isoperimetric problems. The 
idea is to show first the desired inequality for nearly spherical sets following ideas of [22] and 
then use the regularity theory for minimal surfaces to reduce oneself to this situation. To the 
best of our knowledge, it is the first time that this strategy has been implemented to control 
the instability of the ball. Let us also point out that one of the ingredients in our proof is the 
following isoperimetric inequality: 

Lemma 1.7. There exists a constant c > 0 such that for every measurable set E C satisfying 
< al2, there holds 



In recent years, there has been an increasing interest in studying isoperimetric problems with 
densities (see [Ml Emil] for instance). However, most of these authors consider either problems 
where the same density is used for weighting the volume and the perimeter or weights which are 
increasing at infinity. 


In the last part of the paper we come back to the situation g > 1 but consider an infinitely 
stiff trapping potential. That is, we assume that V is equal to zero inside some given open set H 
and is infinite outside. This is somehow the setting which is considered in [421143j . In this case, 
it is easier to work with slightly different parameters. After a new rescaling and some simple 
algebraic manipulations (see Section jl]), the problem can be seen to be equivalent to minimizing 

Jeiv) = e / |Vt/iP + A2|V7/2p + - [ + hi - 1)^ + {K - l)vlvl (8) 

Jn ^ Jn ^ 

with the volume constraint 

/ rjl = ai and / ??i = 02 , 

Jn Jn 

for some A < 1, iL > 1 and ai,a 2 > 0 such that 01 + 02 = |^|- The main difficulty in 
studying the T—convergence of ([ 8 ]) to a sharp limit model is to obtain strong compactness for 
sequences of bounded energy. In the symmetric case A = 1, which corresponds to <7 = 1 in ([1]), 
one can follow the strategy of mm and use a nonlinear sigma model representation |28| to 
rewrite the problem in terms of an amplitude and a phase. In these unknowns, the functional 
takes a form similar to the celebrated Ambrosio-Tortorelli functional |9] from which one can get 
compactness. In the non-symmetric case, the methods of mm do not apply. One has to find 
a different approach. Inspired by recent work on type-I superconductors HZ], we prove instead 
that for every A < 1, the energy ([5]) directly controls a Modica-Mortola [34| type energy where 
r/i and r ]2 are decoupled (see (|M^ L As a consequence, we recover the compactness of sequences 
of bounded energy. Besides extending results of mm to the non-symmetric case, we believe 
that this more direct approach also has an intrinsic interest. In fact, it seems more natural and 
gives a better understanding of how the interaction between the two condensates via the term 
In ^ 1^2 gives to the energy a structure of a double well potential. Using then classical arguments 
from T—convergence such as the slicing method m, we can prove our last main theorem: 

Theorem 1.8. When e —0, T—converges for the strong topology to 




ax,KPiE, Q) ifrji = XE = '^- m and \E\ = oi 
+00 otherwise, 
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where, for a set of finite perimeter E (see |3 [M]/ ). P{E, Q) denotes the relative perimeter of E 
inside and where a\^K > 0 is defined by the one dimensional optimal transition problem 


(^x,K 


= inf 


/R 


Wi\^ + >^W + 


J (?7? + -1)^ + (^ - • 


lim?]! = 0, lim?7i = 1 

— OO H-oo 


At last, in the spirit of what was done in [5S] in the case A = 1 (see also the recent paper 
m), we study the asymptotic behavior of ax^K as K goes to one (mixing) or K goes to infinity 
(strong segregation) and recover good parts of what is expected from the physics literature 

[miiaiiniii], namely 


lim 


<7a,k 


K^i y/K - 1 


2 1 - 
3 1 - A2 ’ 


(^X,K — 0'X,oo 



r\j 


1 aV2 \ 

Kl/2 + ^ 1/4 j ’ 


see Propositions 14.61 and 14.91 

Before closing this introduction, let us point out that most of the results, in particular The¬ 
orem 11.11 Theorem 11.31 and Theorem 11.81 can be easily generalized to arbitrary dimension and 
arbitrary radially symmetric strictly increasing confining potentials. 


The structure of the paper is the following. In Section O we study the case of different 
intracomponent strengths in the presence of a confining potential. We then study in Section [3] 
the crossover case and finally in Section HI we investigate the case of non-equal intracomponent 
strengths in the absence of conhning potential and in a bounded domain. 


Notation 

For X G and r > 0, we denote by Br{x) the ball of radius r centered at x and simply write 
when a; = 0. Given a set E C R^, we let xe be the characteristic function of E. For any integer 
k, we denote by T-& the /c—dimensional Hausdorff measure The letters, c, C denote universal 
constants which can vary from line to line. We also make use of the usual o and O notation. 
The symbols >, < indicate estimates that hold up to a positive constant. For instance, f ff, g 
denotes the existence of a constant C > 0 such that / < Cg. Throughout the paper, with a 
small abuse of language, we call sequence a family (u^) of functions labeled by a continuous 
parameter e G (0, 1]. A subsequence of (u^) is any sequence (ug^.) such that —>• 0 as A: —>• -|-oo. 
For 1 <p < -|-oo and a function /, we denote the norm of / by ||/||p. In the sequel, when it is 
clear from the context, we omit to indicate the integrating measure. In particular, all integrals 
involving boundaries are with respect to and all integrals involving sets are with respect to 
the Lebesgue measure. 


2 The case of different intracomponent strengths 

2.1 The Thomas-Fermi profile 

Let us consider the Thomas-Fermi approximation and prove Theorem 11.11 Let a = (cei, 0 : 2 ) £ 
(0,-|-oo)^. We denote by X the set of pairs of measurable functions (/9i,/02) : R^ R+ and 
by Xa the subset of pairs p £ X satisfying 



aj for j = l,2. 
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For p ^ X we define the energy 

^{p) = [ i;Pl + +Kpip2 + y{pi+P2)- 

7r2 2 ! 2 

As stated in Theoreni ll.il we study the minimization of E in the class Xq,: we prove the existence 
of a unique minimizer for this problem and a stability result. 

We start by studying the minimization problem in the following subsets : given r > 0, let 

= {p^ Xa ■■ supp(/9i) C Br, SUpp(p 2 ) C \ }. 

and 

= 1J^^• (9) 

r>0 



Figure 1: Minimizer of the Thomas-Fermi energy. 


Lemma 2.1. For K > y/g > 1 and r > 0, E admits a unique minimizer in 


Xa’’- Let us set 


r = min(r, ri) with ri = 




and then R? ^ 

TT / 2 2r^ \ TT J 


The unique minimizer is given by 

o,r/ X \{Ri- kn+ */kl < A 

P{ [x] = < 


nr, . I 0 if \x\ < r, ^ ^ 


PYix) = 


^ 0 if\x\>r, ' ' [^{R^-\x\^)+ if\x\>r. 

Optimizing in r > 0, the minimum of E in Xj^ is reached at p^ = p^^o ^ggg Figure\^ with 


ro = ri 


1 + --J-] G(0,n). 

cni y cni J 


In particular, there is a positive gap between p^ and pj at the frontier : 

0-+ = infp? = y/gsnppl > SUPP 2 = cr-- 

BrQ 

The minimal energy is 

Eo = E{p^) = (^{ai + a2f/‘^ + {y/g - 1) a^') ■ 


( 11 ) 


( 12 ) 


(13) 
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Proof. 

Let r > 0. The optimization of E in Xa^ splits into two independent optimization problems. 
From the associated Euler-Lagrange equations, we see that the minimizers have the form given 
by (UnD with Ri, R 2 hxed by the conditions J^2 Pj = aj for j = 1,2. 


Let us now study the variation of r t—)■ = min^^o.r E 

notation t = (r/ri)^, we obtain by direct computation. 


+ E 2 {p 2 ''). Using the 


E^ 


/{*) +1 




with 


fit) 


+ a2rlt 

TT ^ 9 

—rV + a2rn 
o 


for t < 1, 


for t > 1. 


We have / G C^(0,+ 00 ) with f" > 0 in (0,1), f{t) —)■ +00 as t | 0 and f = /'{I) > 0 on 
(1, + 00 ). Therefore, / admits a unique minimiser to = (ro/ri)^ G (0,1). 

Eventually, after some algebraic computations, we get that tq is given by m, we also obtain (fT^ 
and (fT^ . □ 


We are now ready to prove Theorem 11.11 

Proof of Theorem \l.l[ Let p G Xa- We want to show that E{p) > E{p^). Since p^ and E{p^) 
do not depend on K and E{p) is a non decreasing function of K, we only consider the worst 
case K = Then, we write, 

p = p^ + 6 p, with 6 p = { 6 pi, 6 p 2 ). 


The energy expands as 


E{p) = E{p^) + Li{ 6 pi) + L 2 i 6 p 2 ) + Qidp), 


with 

Li{Spi) = f {pi +y^p 2 + \x\^)Spi, 1 . 2 ( 402 ) 


/R2 


(5/0° + V5P? + kP)4o2, 


Qi^p) = l [ ( 4 oi + ^5402)^- 

^ 1 r 2 

This last term is obviously non-negative. Let us study the terms Li{ 6 pi) and 1 . 2 ( 402 )- Let ro, 
l?i, R 2 be as in the dehnition of p^ in Lemma 12.11 We denote by A the annulus \ and 
by U the exterior domain \ Br.^. Using the definition of p\ and P 2 , we have 


Li(49i) = RI [ 6pi + 
J Brn 


^^2 + 

fi- k|2 

IVa ^ 

^ VaJ \ 


4oi + 


'U 


Since fj^2 4oi = 0, we have 4oi = ~ Iauu Using this, we get. 


Li(4oi) = [ ci(|x|)49i, withci(r) = < 


(Rj-r^)-iRl-r^ 


xp4oi- 

for ro < r < 1 ? 2 , 
for r > R 2 . 


AUU 
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r^ — R^ 






We notice that ci is continuous on [ro, +oo) and of the form a + with 6 > 0 on the two 
intervals (ro,ii2) and {R 2 ,+oo). Moreover ci(ro) = 0. Hence ci > 0 in (ro,+oo), ci(ro) = 0, 

c'^(ro) > 0 and ci(r) ’'^3° Now recall that Pi + ^ 0 and since = 0 in H U U, we have 

(5/9 i > 0 in H U U. Hence, Li{6pi) > 0. 

Similarly, 


L 2 {Sp 2 ) 


IBr^UU 


C 2 {\x\) 6 p 2 , 


with C2(r) 


1 r^-Rl 


for r < ro, 
for r > i?2- 


In the unbounded interval {R 2 -,+oo)^ we have C2(r) > 0 with C 2 {R 2 ) 

C2(^) In the interior interval (0, ro), we have C2 > 0, with C2(ro) 

C2 > 0 in [0,ro). Since 5p2 > 0 in Br^ U U, we have -^2(4^2) > 0. 

To sum up, we have establish that for p = p^ + 6p ^ Xa, there holds 

E{p) = E{p^) + H{6p) > E{p^), 

where we have set 

H{ 6 p) = [ ci{\x\)5pi + [ C 2 {\x\) 6 p 2 + 7 ; [ {^Pi + VgSp 2 f- 

JAuU JBroUU ^ Jb? 


= 0, C2(-R2) > 0 and 
= 0, C2(r) < 0. Hence 


Now, for proving the stability estimate ([5]), we have to show that there exists C > 0 such 
that for every function 5p = (dpi, 5 p 2 ) such that p° + dp G there holds 

lldpllii < CH{ 6 p). (14) 

We prove (fTTli in several steps. First let us set i? = i ?2 + 1, \ Br and A = AVP U \ V 

so that BrQ UA[JV = R^. We estimate successively ||dp||i;,i(i/), HdpH^^i^^^ and ||dp||i,i(B^^). 

Step 1. Estimating ||dp||i^i(y). In the set V, we have dpi,dp 2 > 0. Using the fact that ci and C 2 
are positive and increasing on [R,+ 00 ), we deduce 

+ ^ min(ci(fl),C2(Ji)) +<^^(1^1)^ S 

Since the condition p*^ + dp € X^ implies fj ^2 |dpj| < 2aj for j = 1, 2, we get 

lldpllii(^) < Cff(Sp). (15) 

Step 2. Estimating ||dp||j;^i^^^. Notice that dpi > 0 in H. We split the annulus A into three 
subsets A = AiU A 2 U A^ with 

Ai = |x G H : dp2(x) > o| , 

A 2 = |x G A : dp 2 (x) < 0 and dpi(a:) > — 2 y/g 5 p 2 or — y/g 6 p 2 {x) > 2dpi(x)| , 

A 3 = {x G H : dp 2 (x) < 0 and dpi(x)/2 < —y/g 5 p 2 {x) < 2dpi(x)}. 

Step 2.1. In Hi, we have |dpi| + |dp 2 | < dpi + yfg 5 p 2 . Hence, by Cauchy-Schwarz, 

||dp||2,(^^) < |Hi| / {5pr + ^g5p2f < CH{5p). (16) 

J A\ 




Step 2.2. In A 2 , there holds |hpi| + |(5/921 < C\dpi + y/gSp 2 \ and as in the previous step, 

< CH{6p). (17) 

Step.2.3. In we deduce from Spi < —2y/g6p2 and the condition p 2 + (5/92 > 0 that 

Spi < -‘^y/gSp2 < 2y^maxp^ = A. 


Let us note 


m = f 5pi < Oil. 
JAs 


We have 


H{5p) > / ci{\x\)5pi > inf < / ci(|x|)r; : n G [0, A]), v = m>. 


I A3 

Since ci is radial and increasing, the solution of the optimization problem is given by ?;* = 
A1 _Bj,^\_b,,jj where r* > tq is such that = m. Since ci(r) > c'^(ro)(r — tq) for r > vq with 

c'i(ro) > 0 , we deduce 

[ ci(|x|)n*(x) > cm 2 , 

JA 

for some c > 0 depending on A, vq and c'^(ro). This yields 


'9I3 


6 pi ) = m^ <C ci(|x|)c*(x) < CH{5p). 


IA 


Since \5p2\ < in A 3 , we conclude that 

||(5/9||2,(^^^ < CH{6p). 

Gathering (fTHll . tfTTll . lfTHIl . we get 

||(5p||2,(^^ < CH{5p). 

Step 3. Estimating ||( 5 / 9 ||ii(s^^). Proceeding as in Step 2 and exchanging the roles of ci and 
C 2 and 5pi and y/g5p2, we obtain 


(18) 


(19) 


\MIhb.a < CH{5p). 


( 20 ) 

□ 


Eventually, (fT^ . (fT9]l . (f20]l yield the desired estimate (fTTjl . 

Remark 2.2. One cannot hope for a stronger inequality of the form 

\\P - P°lli < E{p) - E{p^) or Up - p^Wl < E{p) - E(p0). 
since one can easily get a contradiction by exchanging p^ and p^ on small balls close to dBr^. 


2.2 Approximation of the Thomas-Fermi limit by the Gross-Pitaievskii func¬ 
tional 

We prove that minimizers of converge to (\/^, V^^) as stated in Theorem 11.21 
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one 


Proof of Theorem \l.SX Let rj’^ be minimizers of Fg, for e G (0,1). Regularizing \/^)) 

can easily construct a competitor % with F^[fi^) — E^je < C. In particular, for minimizers rf^ 
there holds 


e 


< Feiri^) 



C 


so that ([5]) implies 


Wiivif Avlf) - {pi, P2 )\\f < 


Now using that \^/a — y/b\^ < \y/a — Vh\\y/a + Vh\ = |o — 6 |, we have indeed that 





(ll(W)t('7l)")-(p;,P§)lli.)'''" < Ce'l\ 


□ 


3 The crossover case 


We now study what happens when g = \ + for some ^ > 0. For g satisfying ([2|), the energy 
then reads 


F,;{gi,g2) = e [ |V7?ip + |Vr?2|^ 

JR 2 


+ 


^ 2 


^ Vi + lv2 + Kglgl + V{gl+ gl) 


2 Jr 2 


Let us first rewrite the energy in a more convenient way. For this, let g^ be the minimizer of 


GM = e [ \Vg\ 

./R2 


+ 


-f 

£ Jr^ 


g^ + Vg"^ 


under the constraint 



where a = ai + a 2 - It is well known [MCI] that gl converges to 


15={R^-V{x))+ 


where R is such that fj ^2 P = 1- We denote by V the support of p (which is Br when V = \x\‘^). 
We first rewrite the energy in a more convenient form. 

Proposition 3.1. For u = {ui,U 2 ) a pair of non-negative functions, let 


F,iu)=e[ rl2(|Vni|2 + |Vu2|2 ) + - 
Jr 2 £ 


'R 2 


-rf^ (1 - {ul + ul))^ + {K- l)glulul 


then if g = g^u, there holds 


Fe{g) = Ge{g,) + Fe{u) + [ ^g^u^ 

Jr 2 2 

Proof. The proof follows as in |H. We use the Lassoued-Mironescu trick |3D] and write gi = g^ui, 
g 2 = gg,U 2 to get 

|V7?i| 2 + |Vry 2 p = {ul + ul)\Vgf+%V% • V(u? + ul) + gl{\^ui\^ + |Vn 2 p) (21) 
The function 77 ^ solves the Euler-Lagrange equation 

- + \{pl + y{x)Pe) = KPe (22) 
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where Ag is some constant. Multiplying the equation ([2^ by + u^}, integrating and using 
integration by parts we get 

e[ {u‘i + ul)\Vri^\‘^+ ri^Vri^-V{ul + ul) = \ea-- [ + u^)+ Vrij{ul + ul). 

yR2 e Jr2 

On the other hand, multiplying (I22h by fj^ and integrating, we hnd 

X^a = Gs{fj^) + -^ f 

Jr2 

so that (EU leads to 

^X,2 + = ^^£h£(|VuiP + |Vu2p) + ^ { 2 ^^^ ~ 2(^1 +^2)) - VvU'^l+uli^+Gei'ne) 


and therefore 


F, 


{m,m) =Ge{r]^) + f ^rf^uf + sr]l{\Vui\'^ + \Vu2\‘^) 

+ “ (^ 2 ^ ^^ 1^2 + (^1 + ^ 2 )^) + {K — l)rj'^ufu 2 ^ 

= Ge{^e)+ [ + S f ?72(|Vmi|2 + |VU2|^) 

./r 2 ^ Jr 2 


1 

+ - 


£ JR2 V2 


(1 - (uj + ul)f + {K - l)rf^u\ul 




= Ge{r]e)+Fe{ui,U 2 )+ -'n^U 2 

J-R? 


2 ■'£ 2 ) 


□ 


aKF{E) if Ui = xe,U2 = Xe<^ and V{E) = ai, 
+00 otherwise. 


which completes the proof. 

For E a set of locally hnite perimeter in P (see El EH), let 

F{E) = [ -fG and V(F;) = [ P- 

JdE Je 

It is proved in [2S1 Th. 1.1] (see also IH) that for all p < 00, Pg L^ —r converges to the functional 

Qo{u) = 

where ax > 0 is dehned by the one dimensional optimal transition problem 

ax = inf jp'ip + |? 7^|2 + 1 (ry2 + 7^2 _ ^^2 _ i)y^2^2 . ^ ^ 

Since | f is a continuous perturbation of E^{ui,U 2 )j we immediately obtain the following 
result. 

Theorem 3.2. For every p < 00, the functional Fi.{ui,U 2 ) + 1/ — T converges to 

{axF{E) + ^j if ui = xe,U 2 = XE‘^ andV{E) = ai, 


Qi{u) = 


\ 


+00 


otherwise. 
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Up to dividing by ax and modifying we can assume that ctk = 1 - If ui = xe, U 2 = Xe<^, 
we will, by a slight abuse of notation denote G^{u) by Q^{E). We now want to study the 
minimizers of (whose existence follows from the Direct Method). As in [2T], by making a 
spherical symmetrization, we can restrict the analysis to spherically symmetric sets. For a given 
set E and a given half line £ starting from zero, such symmetrization is defined by replacing for 
every r > 0 the spherical slice E n dBr by the spherical cap K{E,r) centered in n dBr and 
such that 'H^{K{E,r)) ='H^{EndBr). 

Proposition 3.3. For every minimizer E^ ofQ^, dE^ DV (recall that T) is the support ofp) is 
a C°° hypersurface. Moreover, for every half line i starting from zero, there exists a minimizer 
E^ such that for every r > 0, dBr n E^ is an arc of circle centered in £. 

Proof. Let E^ be a minimizer of Q^. The regularity of dE^ is a consequence of the regularity 
theory for minimal surfaces. Indeed, since p is locally bounded away from zero in B, any 
minimizer of is locally a quasi-minimizer of the perimeter (see [211 Th. 3.2] or |31[ IT9] 
for instance). From this, one can infer (7^’“ regularity of E^. Since p is smooth in B further 
regularity follows. The symmetry follows as in [211 Th. 3.2] since spherical symmetrization 
reduces E and leaves invariant. □ 

When ^ is small, the perimeter term is dominant. In this case, we are in a situation similar 
to the one studied in mm- 

Proposition 3.4. There exists oq G [0,1/2) such that for every ai G (ao,a/2 —oq]; there exists 
^o(ai) such that for every ^ < ^o(ai), the minimizer of is not radially symmetric. 

Proof. It is proved in [1] that for every such ai, the minimizer of E under volume constraint is 
not radially symmetric and thus there exists with 

E{Eaf) < inf {.F(F) : F radially symmetric } . 

V(F)=q:i 


but since 

G^{Ea,) = EiEr,,,) + ^ [ f 

it is clear that 

G^{Eaj^) < inf {E{E) : E radially symmetric } < inf{^^(F) : F radially symmetric } 
V(F)=oi 

for ^ small enough. □ 

We now study the situation of large Our main result is a rigidity result stating that for 
large (but not infinite) the unique minimizer is the centered ball. 

Theorem 3.5. For every ai G (0,a), there exists .^i(q;) such that for every ^ > ^i(q;), the 
unique minimizer of is the centered ball B such that V{B) = ai. 

In the rest of the section, oi is fixed. In order to ease notation, we assume that the unit 
ball is such that V{B) = ai (the general case follows by dilation). Theorem 13.51 follows from a 
combination of two results. The first one is a stability result for the volume term 

Proposition 3.6. There exists c > 0 such that for every set E with V{E) = V{B), 
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Proof. Let E be such that V{E) = V{B) then 


Let 


and 5 be such that 


[ [ 75(1 - Ixp) + [ p(\x\^ - 1) 

Je<= Jb^ JE^^nB JB^nE 

L 

j> = v. 


V= I p 

lE^nB 


>B\Bi_s 


Letting Fs = B\Bi_s, since f^E<=r,B)\Fs P = /fas" P' 


[ pa-\x\^)- [ pa-\xf)> ( 

JE^nB JFs \ 


inf (1 — |xp) — sup (1 — |xp) 


(F=ns)\F^ 


FAF" 


IFs\E‘^ 


P>0 


and therefore, Fs minimizes p(l — |xp) among sets G d B with JqP = V so that 




P (1 



xp) > cV"^ 


(24) 


Similarly, letting Eg = n B^ with Jp_p = V = Jpcf^pp, then F^ minimizes Jq p{\x\‘^ — 1) 
among G <Z B^ with f(^P = V and thus 


'B'^nE 


p(lx|^ — 1) > p{\x\^ “ 1) > CV 


Together with (IMll . this gives (1^ . 


□ 


The second is an estimate on the possible instability of the ball for E. 

Proposition 3.7. There exist e > 0 and G > 0 such that for every set E with V{E) = V{B) 
andJp^pP<e ^ 

F{E)-F{B) >-g([ p) . (25) 

\Jeab J 

Since the proof of Proposition 13.71 is long and involved, we postpone it. Let us show first 
how Proposition 13.61 and Proposition 13.71 yield together Theorem 13.51 Let E^ be a minimizer of 
then using B as competitor, we obtain thanks to (j23p and (|25p . 



1 i[g^{E^) - F{E^)] - [g^{B) - F{B)]) 


< 1{F{B)-F{E^)) < 


C 

cf, 



This implies 


(L 


EjAB P 


= 0 


for G jc and concludes the proof of Theorem 13.51 


Before going into the proof of Proposition 13.71 let us comment a bit on the statement and 
give the strategy for proving it. As explained below, in general, we do not expect the ball to be 
a local minimizer of F (and in particular, we cannot expect (I25p to be true with a plus sign on 


13 


the right-hand side). However, (I25p shows that in some sense, the Hessian of F at the ball is 
bounded from below. The proof of (I25p follows the strategy of m for proving the quantitative 
isoperimetric inequality (see also mm)- Inequality (j25p is first shown for nearly spherical sets, 
borrowing ideas from [22]. The proof is then finished by arguing by contradiction, construct¬ 
ing a sequence {En) converging to the ball and contradicting the inequality. Using a Selection 
Principle and regularity theory for minimal surfaces, it is possible to replace (En) by a better 
sequence (En), still contradicting (|25p but converging in a much stronger way to the ball. Since 
for n large enough En are nearly spherical, we reach a contradiction and ProDosition l3.7l is proved. 


We thus start by proving (I25p for nearly spherical sets (in the sense of Fuglede [22] ). 

Proposition 3.8. There exist eo > o.nd Cq > 0 such that ife < Eq, every set dE = {(l-|-u(x))x : 
X G dB} with ||u||Lip < e and V{E) = V{B) satisfies 


F{E) - F{B) > -Co 


'eab 


Proof. Recall that p{x) = {P? — |xp)+ with R> 1. The condition p = f^p can be written 
as 

JdB ^ 4 


which leads to 


L 


13-R2 


u = 


I OB 2 R2 _ 1 
We can now compute the energy 


/ 

JdB 


u -b O llu 


iLRas )) ■ 


(26) 


(R2 _ 




2 u + u^\ 
R2 _ 1 


(1 + u) 1 + 


|Vu|- 

(1 -b ufi 


2 \ 1/2 


- 1 


'dB 


U 1 - 


R 2 _ 


l) 2 (R2 _ 1)2'“^ + IfW'^'^Wl'^idB) + O (11111112(96) + ll^ll|li2(9B)) 


Using 


this turns into 


(R2 _ 1)3/2 “ -^(^)) - 2 + ll^“lli2(9S) + "(ll^lliRaB) + ll^'“lli2(9S)) 


> o Il^ii|li2(a6) C'l|ii|li2(96)- 


(27) 


We now claim that for every <5 > 0, there exists > 0 such that 

2 


[ < 6 [ |Vu|2 + A5 ( 

JdB JdB 


\U\ 


dB 


(28) 


Indeed, denoting by u = ^ u and using Sobolev embedding and Young inequality, we 
compute 


/ = / (u-u)u-b27r(u)^ < ||u-u||x,oo( 96 )||ii||l 1 ( 9 b)+ (2vr) 

JdB JdB 

~ II^1i||l2(96)||u||61(96) + I|ii|Il1(9B) 

< <5||Vt(||^2(96) + ^<5l|ll|lii(9B)- 


-1|L.I|2 


\B{dB) 
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This proves (I28I1 . Using (f28]l in (f27)l . we obtain, 

T{E)-T{B) > -c( [ [ 

\JdB 

Now, since for nearly spherical sets, there holds 

/ p ~ \EAB\ ~ / |ri|, 

JeAB JdB 

we can finally conclude the proof. 

Remark 3.9. In the proof of Proposition 13.81 if we assume that E is centered, that is 

/ 'p{x)x = 0 , 

Je 

we could do as in [221 [IS] and decompose u in Fourier series on dB to get 

M = ^ UkYk 


□ 


with Yq = 1 and Yi = x ■ v for an appropriately chosen v £ dB. So that 


I 




1 


1 




1 


We find thanks to 


uo= u = 0{\\u 
JdB 


ui = f uYi = O (\\u 
JdB 


From this, letting = k‘^ being the k—ih. eigenvalue of the Laplacian on dB, 

^2^0_L^^2^ r — R2^2 + R2) 

k>2 

Since > 4 for k > 2, and since 


X ■ E = Q. 


'dB 


f R^{2 + R^) 2 /■ i2 V- 2/ + 

Jm ~ (R^ - If “ )+»{ll“ll2) 


4_UP±^>o 

(R2 _ 1)2 


1 (5_V13)V2 

OT R > ±{5 + Vu)^/^ 
v3 v3 


we expect the ball to be unstable for R G ((5 — \/l3)^^^/\/3, (5 + \/l3)^/^/\/3) and stable 
otherwise. 


In order to go further, we need an isoperimetric inequality which is, we believe, of independent 
interest. 

Lemma 3.10. There exists c = c{a) > 0 such that for every set E with V{E) < a/2, we have 

F{E) > cV{Ef/^. (29) 

Proof. Up to a dilation, we assume here that R = 1 that is ^(x) = (1 — l^cp)^. Since the estimate 
is sublinear, it is moreover enough to prove it for connected sets. Let us first show that we can 
assume that E C [Bii^Y- Indeed, if this is not the case then there exists x G dE n Then, 

there are two possibilities. Either E C R 3/4 and the classical isoperimetric inequality already 
gives ([29]l or there is y G dE n R 3 / 4 . Since dE is connected it contains a path from x G R 1/2 
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to y e s-iid J^{E) is bounded from below by — s^) = 29/192. Again, this gives (f2U|l . 

As in Proposition 13.31 we can make a spherical rearrangement and assume further that E is 
spherically symmetric. 


We now transform our problem in order to work on a periodic strip. Let be the unit torus 
and consider the diffeomorphism (j) : Bi\Bi /2 —>• 5"^ x ( 0 , 1 / 2 ) given by 


(/(6l,r) 



For E C X (0,1/2), we let 

n l/2 

XF{x,y)y dydx 

then, since for E C Bi\Bi/ 2 , 


V{E) r^Vi{ct){E)) and 


and ^i(F) = [ y3/2 
JdF 


E{E) EiicPiE)), 


we are left to prove 

Vi(F) < yFcS^x (0,1/2). (30) 

Notice also that since the sets E C Bi\Bi /2 started with were spherically symmetric, we 
can further assume that for every y G (0,1/2), F r]{y = y} is a segment centered in 0 x {y} for 
instance. We then make a convexification step. We define F as the smallest set which contains 
F and which is y-convex, that is for very x £ the set {y > 0 : {x, y) G F} is a segment (see 
Figure [2|). We have 

y A 



Figure 2: Convexification step in the y-direction. 

Vi(F) > Vi(F) and Fi(F) < Ei{F). 

Indeed, the first inequality comes from F C F and the second from the fact that for every x G S^, 
the shortest (weighted) path between the points (x,yi) and (x,y 2 ) is the straight segment. It 
is therefore enough to prove (1301) for the sets F. After these two symmetrizations, the set F 
is contained between two graphs yi : (—T, T] —)• (0,1/2) and y 2 : (—T, T] —)• (0,1/2) for some 
0 < T < 1/2 i.e. 

F = {{x,y) : -T <x <T, yi{x) <y< y 2 {x)}. 

Moreover, yi and y 2 are even, yi is non-decreasing and y 2 is non-increasing in [0, T]. Using this 
parameterization, we have 

Vi(F) = l r [y2 - y2] and Fi(F) = y^^^^l + \Dy,\^ + yf' y^TTJ^, 

^ . -T . -T . -T 
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where for a function y of locally bounded variation, 


Ti{y)= f + \Dy\'^ ■■= [ y^/‘^^/lTW?+'^\y{x'^f/‘^-y{x y^/‘^d\D^y\, 

x&jy 


see P Th. 5.54], Notice that we can assume that yi(±T) = y 2 {±T). Indeed, if T = 1/2, then 
J^i(F) > ^ 2 ( 1 / 2 )^^^. In this case, we can add a vertical cut between ^ 2 ( 1 / 2 ) and yi(l/2). The 
additional contribution to the weighted perimeter is of order of ^ 2 ( 1 / 2 )^'^^ — yi(l/2)®/^ which is 
controlled by J-i{F). Eventually, since Fi is the relaxed functional of its restriction to paths 
with respect to convergence [SI Th. 5.54], we can further assume that yi and y 2 are smooth. 
After, these symmetrization and regularization steps, we can write, 

Fi{F) ~ r yl^^ + r + ^ 2 ( 0 )^/^ - yi{0f/^. (31) 

J-T J-T 

We then consider two cases. First, if yi(0) <C ^ 2 ( 0 ), then (l3T|) yields 


from which we deduce 

Vi(T) > r yl < 2/2(0)'/" r vT < T^{Ffl^T^{F) < T^{Ffl\ 

J-T J-T 

Then, in the case yi(0) ~ 2 / 2 ( 0 ), by convexity 2 / 2 ( 0 )^/" — 2 / 1 ( 0 )'’/" > |yi( 0 )'^/"(y 2 ( 0 ) - 2 / 1 ( 0 )) so 
that (fHTT) becomes 

-^i(^) > Tyl/\^) + 2/2(0)3/"(y2(0) - yi(0)). 

We infer that 

Vi(F) = y^(2/2-2 /i)(2/2+ 2/1) < 4 ry 2 ( 0 )(y 2 ( 0 ) - yi( 0 )) 

< 4Ty2(0)y2(0)'/®(y2(0)-yi(0))'/5 = 4ry2(0)''/" (2/2(0)''/"(y2(0) - yi(0))) 

< Fi{F)Fi{Ff/^ < Fi{Ff/\ 


which concludes the proof of (130]) . 

□ 

Remark 3.11. The exponent 5/6 in (1291) can be easily seen to be optimal by considering as 
competitor a small ball touching the boundary of B. 

We can now prove the following e-regularity result: 

Proposition 3.12. Let A > 0. Then, there exists ei > 0 such that if E is a A-minimizer of F, 
i.e. for every G, 

F{E) < F{G) +A [ -p 

Jeag 

and if < £1 then E is nearly spherical i.e. dE = {(1 + u{x))x : x G dB} and 

||it||ci-“ ^ ^0 (where Eq is the one defined in ProDosition \d. 8\} . 
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Proof. Fix 5 > 0 then inside -B 1+5 by classical e-regularity results for quasi-minimizers of the 
perimeter (see [T9l Th. 6.1] or [SU Th. II.6.3] for instance as well as [HI Lem. 3.6]), if 
I{EnBi+s)^B'P — then dECiBi^s is a small (7^’“ perturbation of dB in particular, if < 

ei, dE n i?i +5 C B,.s and E can be written as E = Ei U E 2 where Ei is nearly spherical, 
E 2 C By testing the A—minimality of E against Ei, we find by (|29p 

cViE2)^^^ < F{E2) < KV{E2) 

and thus if V(£' 2 ) / 0 , 

(^) <V{E2 )<ei 

which is absurd for ei small enough. □ 

We will also need the following simple lemma which is a weak version of (|25p . 

Lemma 3.13. There exists Ai > 0 such that for every set E, 


F{E) - F{B) > -Ai 


lEAB 


P- 


(32) 


Proof. Let n be a vector field with suppn C Br, |n| < 1, n = x on dB and jjdivnjjoo < C and 
let 

w = • V -|- p^/^divn, 

so that div (p^/^n) = ptc. Then, for every set E, denoting by the outward normal to E, 


F{E) - F{B) > [ 
JdE 


JdB 


/ pw — / pw > —lire] 

Je Jb 


00 P- 

JEAB 


□ 


We can finally prove (|25p . 

Proof. We argue by contradiction. Assume that there exists a sequence of measurable sets (F'n) 
^ith Je^abP^^ and 

F{En) + 2Co( [ p) <F{B), (33) 

\JEnAB J 

where Cq is the constant given by Proposition 13.81 Let En = as P and for Ai > 0 given by 
Lemma l3.131 and A 2 > 0, let Fn be a minimizer of 


F{F) + 2Ai 



1 2 


I FAB 


+ £r, 


1/2 

+ A 2 lV(F)-V(B)l. 


(34) 


Step 1. We claim that if A 2 is large enough, V(F„) = V[B) for all re. 

To prove this, we follow the approach of [20] (see also [23] for another approach). Assume 
by contradiction that the claim does not hold, then there exist sequences of positive numbers 
—>• 0 and A 2 ,a: —^ + 00 , and a sequence of measurable sets Gk C such that Gk minimizes 


F{G) + 2 Ai 



n 2 


/GAB 


1/2 


+ 


+ A2 ,felV(G)-V(R)l 


(35) 
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and for instance V(Gfc) < V{B) (the other case is similar). In order to get a contradiction, we 
build a new sequence Gk such that V{Gk) = V{B) and 

\ 1/2 


BiGk) + 2Ai 


^k- L P 

JGk^B . 


+ 


< 


\v{Gk)-vm- 


The construction is a bit delicate since we want all the constants to be uniform in If this 
were not the case, one could have simply used Almgren’s construction (see [31]). 

First, testing the energy (13511 with B, we find that F{Gk) < G and therefore, up to extraction, 
Gk converges in to some Goo- Moreover, since \V{Gk) — V(i?)| < Goo satisfies 

V(Goo) = V{B). The set Goo, minimizes 


.F(G) + 2Ai / p 

Jgab 


under the constraint V(G) = V{B) (this can be seen for instance by a T—convergence argument). 
Therefore, by (1521) . Goo = B. Notice that arguing as in Proposition 13.31 we get that the sets 
Gfc n 2? are G^’". Now, let us fix d > 0 and r <C 1, let us choose xq G dB and let us set 
xo = (1 + ^)xo and G = (1/2) since Bj./ 2 ixo) CiB = 0, there hold 


L 


p < 6 and / p > G (36) 

f Gl^r\B^/2(xo) JGkClBrixo) 

for k large enough. Let 0<(7fc<l/2bea sequence to be fixed later and consider the bilipschitz 
maps: 

(1 — 3ak)x if |x| < 

x + ak{l-^)x if § < |x| < r, 

^x if |x| > r. 


^k{xo + x)-xo= < 


and let Gk = ^k{Ek)- 
Step 1.1. We first prove that 


F{Gk)-B{Gk)>-GakHGk)- 


(37) 


Following the notation of |20|, we let for x G dGk, Tk^x{x) = V<I>fc(x) o r for r G -nk.x (where 
TTk^x is the tangent space to dGk at x) and 


JiTk,x = \J^^^{Tl^x°Tk,x) 
be the one Jacobian of Tk^x so that 

B{Gk) = [ ■f/\^^^x))j^Tk,xd'H\ 

JdGk 

In particular, it is proved in [20] that JiTk^x < 1 in Bj./ 2 {xo) and 

diTk^x ^ 1 T ^^k 

in Gr{xo) = Br{xo)\Bx/ 2 {xq). We can now decompose, 

F{Gk) - HGk) = [ - [ // 2 ($^(^)) 

JdGkr\Br{xo) JdGkr\Br{xo) 


L 


dGknCrixo) 


fG- 


/aGfcnaixo) 


+ 


L 


9GfcnBr/2(a:o) 


fG- 


'9GfcnSj./2(a;o) 
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But since in Br{xo), \p^^‘^{x) — < Cp^^‘^{x)ak, 


IdGknCrixo) 


Jd 


(1 - Jirfc,,)p3/2 

aGknCr{.xo) 

+ JiTn,i{f/\x) - f/\^k{x)))dH^ 
[3/2 


> -C(Tk / P' 

JdGknCr{xo) 


Similarly, 


'dGknB^/2ixo) 


py‘^-py\,^,(x))JiTk,,d'H^= [ 

Jd 


(1 - JiTfc 

9GfcnB,/2(xo) 

> -Cak f 

JdGkr\B^/ 2 {xo) 


fl\ 


Hence, pTll follows. 

Step 1.2. We now prove that for some k > 0, 

V{Gk) - V{Gk) > KUk. 

If J^k denotes the Jacobian of d>fc, it is shown in [20) that 

J^k < 1 + ScJfc 

and that in C'r(xo), 


J^k > 1 + CCJfc 


(38) 

(39) 

(40) 


for some c > 0. We can decompose 


'Gfc 


P = 


'GkCiBr (xq) 


p{^k{x))J^k{x) - p{x) 


JGknGr{xo) 

Thanks to (l39]) and (|40l) . there holds 


p{^k{x))J^k{x) - p{x) + / p{^k{x))J^k{x) - p{x). 

J Gfcni?j./2(a:o) 


L 


GfcnC*r (3I0) 


p(d>fc(x))Jd>fc(x) - p{x) = 


' GfePlOr (iTo) 

> cak 


(p(d>fc(x)) - p{x))J^k{x) + {J^k - ^)p 


/ 

JGtnGrlajn') Jg, 


p -2 j |p(4>fc(x)) - p(x))| 

GfenGr(a:o) JG^nGrixo) 

-p{^k{x)) -'p{x) 


' GfePlCr (iTo) 


P\C(Tk-2 


p{x) 


Since p{x) = {B? — |xp)+, we see from the definition of d>fc that in Gr{xo), 


Pi'^kix)) - p{x) 


p{x) 


< GakV 


hence we can choose r small enough so that 2Gakr < C(Tfc/2. So that 

/ p{^k{x))J^kix)-p{x) > Kiak p, (41) 

J Gfen(7r(lTo) J GkC\Cr{xQ) 
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for some ki > 0. Similarly, thanks to (1401) . 


/ p{^kix))J^kix)-p{x)= (p(4>fc(x)) - p(x))J$fc(x) + (J4>fc - 1)/) 

'GfcnB^/2(2;o) Gkr\B^/2(xo) 


> -K2Crfc / p, 

J Gk^B^i 2 


for some k ,2 > 0. Combining this with (j4ip and (I36p gives 

so that pHp holds if 5 is small enough. 


V{Gk) - V{Gk) > o-fc Ki / P- K2 I p] > (kiC - K2S)ak, 

y JGf^nBrixo) JGkr\B^/2{xo) 


Step L5. Since V{Gk) < V{B) and V{Gk) and the continuity of the map 

fXk h’(Gfc), show that we can find (jfc —)• 0 such that V{Gk) = V{B). From this we get a 
contradiction. Indeed, by minimality of Gk, 


T{Gk) + 2Ai 


UGkAB 


p-^k 


\ 1/2 

+ e'J +A2_,|v(Gfc)-V(S)| 


< B{Gk) + 2Ai 




n2 \V2 
"k 


p — e'k + 


< B{Gk) + Gak + 2Ai 


UGfcAS 


P-4 


2 \ V2 

+ 1 


where in the last line we used (l37P . Since the function x {{x — + 4)^'^^ i® 1—Lipschitz, 

we obtain, recalling 


^ 2 ,k<Xk ^ (Xk + 

Arguing as for pHp . we can prove that 


J~ I ^ 

JGkAB JGkAB 


L p-f p 

JGkAB JGkAB 


4 ^ki 


from which we obtain 

This is not possible since A 2 ^k —^ +oo. 


A2,k(Xk < fXk. 


Step 2. Going back to minimizers of (|Mp . we have that Fn converge to some F^o 
minimizing 

F{F) + 2A, [ p + A 2 \V{F)-V{B)\ 

Jfab 

that is by (1321) . F^o = B. The Fn are A-minimizers of F. Indeed, for every E, 


F{Fn) < F{E) + 2Ai 


/ ^ 

+ Aa 

f,-f. 

JeAB JEnAB 


JFn J E 


< F{E) + (2Ai + Aa) 


lEAFr, 
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Therefore, by Proposition 13.121 and Proposition 13.81 for n large enough, 


F{Fn) - F{B) > -Co 



(42) 


Step 3. Let 7 ^ = p. Then, by minimality of there holds (using En as a competitor 

and recalling (f33|) i. 

E{Fn) + 2Ai ((7„ - En)^ + < F{En) + 2Ai4/2 < F{B) + 2KieY^ - 2Co4 

which combined with (j42l) gives 

2Al ((7n — Sn)^ + £n) ^ ~ C'oTn — + ^Aie^^. 

From this we obtain 

((7n — ^ < ^" 07 ^ + 2 Aiey^. 

Dividing by 2Ai and taking the square of both sides, we get 

{in - + En < 7^7n + 


4A2 Ai 


subtracting and dividing by 7 ^ we deduce 


2 --(2 

1- — 1 <7A7n + 2w^^^ 

7ny 4A2 Ai 


^n\ ^ tyg ^ 2 I 0^*0 ^1/2 


hence 1 — ^ 0. Going back to (|43l) . we obtain 


2A,4-'" 


I — - 1 I + 1 


nl/2 


-1 <-Goe 2 2 _ 




(43) 


from which we get a contradiction since for n large enough the left-hand side is positive while 
the right-hand side is negative. This ends the proof of Proposition 13.71 □ 


4 The case without confining potential 

4.1 The functional 

We finally go back to the situation where g > 1 but where there is no trapping potential. We 
thus consider for a fixed bounded open set D C 

Feiv) = £ [ |V?7ip -f |Vr/2|^ + -[ lit + f72 + ^7i72 
Jn e Jo ^ 2 


with the constraints 


Letting 



7 = (ai + “25^^^), 


(44) 
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it can be easily seen that minimizers of the Thomas-Fermi energy 

1-2 , 9-2 


under the mass constraint 


E{p) — 2^1 + 2^2 + KpiP2 


Pi = Oil 


7 


are given by 

for any set E G ft with |-E| = . The minimal energy is then 


Eo = E{p) = 


7 


This motivates the change of variables 


= (m) 

yielding 


- 1/2 


m, m = 


7 


\n\gy^ 


- 1/2 


m, X 


2\n\ 


mj 


I, J,® = 


= e L I Wit + 9 ‘'"^1 Viftt +1 / 1 (ijt + m‘ -1)’ + ® -1)91 
Jn ^ Jn ^ 

where 1 < K = under the mass constraint 


/ 


~2 ~ 
Pi = Oli 


where ai = ai and a 2 = 9 ^^“^012 so that ai + 02 = 7 = |n|. Forgetting the tildas and letting 
1 > A2 = g we finally obtain that the original minimization problem is equivalent to 
minimizing 

Je{p) = £ [ iVr/ip + A2|VT/2p + - f + + iK - l)pipl 

Jn ^ Jn ^ 

under the volume constraint 


4.2 The one dimensional transition problem 

Let us introduce, the following energy defined for ?/ = (71,1/2) S 

Ex,k{p) = [ \Pi\^ + ><^\V 2 \^ + Wk{pi,P 2 ), 

Jn 

where for s, t > 0, we introduced the potential 

WK{s,t) = ^{1 - + {K - . 

We consider the minimization problem, 

ax,K = inf\Ex,K{p) ■ 7 = ( 71 , 72 ) G IF/c (R,R+), lim 7 i = 0, lim 7 i = 11 . (45) 

I —00 +00 I 

Let us show that problem (|45p admits a minimizer. The result also follows from [a Th. 2.1] 
(see also [7]) with a different proof. Let us point out that uniqueness of the optimal profile has 
been recently shown |40j . 
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Proposition 4.1. There exist minimizing pairs to crx x- Every such minimizing pair is smooth. 
Moreover, the following equipartition of energy holds: 

Proof. We establish the existence of a minimizing pair 77 = (7/1,772) by the Direct Method of 
the calculus of variations. The required compactness and semi-continuity result is stated in 
Lemma Ml The smoothness of 77 is the consequence of the Euler-Lagrange equations 

( - 2 %' ) = 0, in R. 

For the equipartition of energy, we take the dot product of the Euler-Lagrange equations with 
{h'i^v'2)^■ Integrating, we see that the quantity —It/'j^P — A^jr/^p -|- ILa'( 77 i, 7/2) does not depend 
on X. Using the conditions at inhnity, we conclude that [t/'^P -|- A^It/^P = 11^(7/1,772) on R. □ 

Lemma 4.2. Let {r]^)k>o C iy,]^’^^(R, R+) be such that there exists Cq > 0 such that for k > 0, 
hminf77j^ <1/4, hmsup77j^ >3/4 and Ex < Cq. 

OQ 3;_)._|_oo 

Then, there exist 77 = (771,772) G IE/,p(R, R^), a subsequence (still denoted by and a 

sequence (zk) C R such that 

r]^{- — Zk) —>■ 77 uniformly on any bounded subset o/R. 


Moreover, 


and 


lim 77 = (0,1), lim 77 = (1,0), 

X—>■—CXD X^-I-OO 


Ex,K{‘n) < lim inf Ea,1^(77^). 

fc^ + OO 


( 46 ) 


To prepare for the proof, we start by noticing that (similarly to the Ginzburg-Landau func¬ 
tional for Type-I superconductors m) the Gross-Pitaievskii energy controls a classical double 
well potential. This is the key ingredient in order to get compactness for sequences of bounded 
energy. Let us introduce the relaxed potential. 


wk{s) 


inf Wk{s, t) 
tsR 


^(1 -s2)2 - 1(1 -iLs2)2 if0<S<E:-V2^ 

l(l-s2)2 ifs>iL-l/2. 

< 2 


( 47 ) 


The function wk is a standard double-well potential (see Figure [ 3 |). In particular, since there 
holds wk{s) = WK{t) < WK{s,t), we have for 77 = (771,772) G (R, R+) 


/R 


1 771P+ 7771^(771) < Ex,k{v) and 


A 17 ? 2 p+ 7771^(7/2) < Ex,k{v)- ( 48 ) 


7R 

Proof of Lemma Let us consider a sequence 77^ = (77^, 7/2) and Cq > 0 satisfying the hy¬ 
potheses of the lemma. Let us hx A: > 0 and let us consider y < z such that either, T]((y) = 1 / 4 , 
77 ^(^;) = 3/4 or 77^(77) = 3 / 4 , r/^(z) = 1 / 4 . Using the Modica-Mortola trick, we see that 


/•2 /- /■3/4 _ 

I (771 )'p -h 777x7/1 > 2 / 1 / 777x7/1 I (7/t)' I > 2 / x/wk{s) ds = 6. 
Jv '' 4 i /4 
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Figure 3: Relaxed potential wk{s) and limit potential •w{s) 


{1 — (dashed line). 


Since <5 > 0, taking into account (H8|) and the bound E\^k{v^) < Cq, we deduce that there exists 
an odd integer G [1, Cq/5\ and sequence of intervals 

^0 = (-00,4), 4 = (4,4), ••• , 4,-1 = (4,-1,4j, 4, = (a:n,,+oo), 


such that r]^{Xj) = 1/2 for j = 1, • • • ,and there exists Uq,- " 


, 4, with yj G such that 


j ViiUj) = 1/4 and r/f < 3/4 in Ij if j is even, 

= 3/4 and rj^ > 1/4 in Ij if j is odd. 

See the example of Figure [H Now, up to extraction, we assume that Uk = n does not depend 



Figure 4: Example of construction of the intervals Ij. 

on k. Now, if the sequence {x^ — x\)k is not bounded, we can extract a subsequence such that 
(^2 — x\) —)■ oo. Repeating the process, we can assume that the following property holds true: 
there exist R > 0 and a partition of {1,.., n} = {ao,... , Iq} U {oi,..., bi} U • • • U {am, ■ ■ ■ , bm} 

with qq = 1 , bm = n and Oj+i = + 1 for j = 0 ,..., m — 1, such that, 

V/G{0,---,m} 0<4-4i < ^, V/G {O,--- ,m-1} 4^+1 “ 4, +oo. 

Since oq = 1 and bm = n are odd there exists Z* G {0, • • • , m} such that and bi^ are odd. 

We set Zk = x^^ and yf = rif{- — Zk). By construction, there exists a sequence (Rk) C [R, +oo) 

with Rfc —>• +00 such that for k > 0 , 

4 <3/4 in[-Rfc,0], 4(0) = 1/2 and 4 > 1/4 m[R,RA:]. (49) 

Now, from the energy bound Ex^k (4) < Cq, we see that (4) is bounded in l^oc (R). Up 
to extraction, there exists rji G W^oc (R) such that r]i{- — Zk) = Vi ^ Vi locally uniformly. 
Moreover, from (1491) . we have 

limsup??! < 3/4, 771(0) = 1/2 and liminf 771 > 1/4. (50) 

-CXD +°0 
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Similarly, up to extraction, there exists 772 G such that r}2{- — Zk) —)• 772 locally uniformly. 

By lower semi-continuity of the Dirichlet energy, we have 

Ex^xivi^m) Ex,Ki'ni,'n2)- 

k^+oo 

To end the proof of Lemma [ 4 .2[ we have to establish that rj = (771,772) satisfies the conditions at 
infinity (HHIl . Since is hnite 77 admits limits 77^ at ±00. From the bound WK{'r]i-,V2) < 

C, we get 77“, 77+ G {( 0 , 1 ), ( 1 , 0 )}. Eventually (I^Ujl implies 77“ = ( 0 , 1 ), 77+ = ( 1 , 0 ), that 
is dm). □ 


4.3 The F—convergence result 

In this section we stndy the F—convergence of Jg as e goes to zero. 
Theorem 4.3. When e ^ 0, Jg T—converges for the strong topology to 




ax,KP{E, Q) i/ 771 = XE = 1 - 7?2 and \E\ = ai 
-boo otherwise. 


Proof. Since Jg(min(77i,l -|- 5 ),min(772 ,1 -|- 6)) < Jeip) for all 5 > 0 and since the bound 
In £ + (^2)^ — 1)^ < C* implies that if min(77|, 1 -|- ( 5 ) (strongly) converges to some 77 then 

also 77g converges to the same 77, we can always assume that sequences which are bounded in 
energy are bounded in L°°. 

Now, we use again that the Gross-Pitaievskii potential controls the double well-potential wjc 
(see dm)- We have, 

/ WK(m) + / yjKim) < 2 / ^^(771,772). ( 5 i) 

J ^ J Q, J Q, 

Thanks to m and the usual Modica-Mortola argument, we then have that from every 
sequence (771,77!) of bounded energy, we can extract a subsequence converging strongly in to 
some pair (771,772). Moreover, from the bound on the energy, we get that 771 (x), 772(3:) G { 0 , 1 } and 
7?i (37)772(x) = 0 for almost every x G R. From the volume constraint and the strong convergence, 
we also deduce that rjf = ai. The lower bound inequality is then a standard application of 
the slicing technique (see [13 H ES] for instance). The upper bound is also standard. By 
approximation it is enough doing the construction for a smooth set E. Let 5 > 0 be hxed then 
we can find T > 0 and (77^,772) with 7jf{—T) = = 0 and r]f{T) = p^i—T) = 1 with 

J + "^^1(^2)^^ + ^^-(77^,772) < ax,K + d 


Let he the signed distance to dE. 
We then let 


rf(x) = r]^ 


dEix) 


Using the coarea formula it can be seen that 77^ converges strongly to {xe, 1 — Xe) and that 


limsup Jg(77'^) < {ax^K + d)P{E,Q). 
£—^0 


Letting finally 





we have that by definition satisfies the mass constraint. Moreover, using that \\r]f\\2 = cti+ 0 {e) 
(and therefore, ry? = ryj(l + 0(e))) we have 

limsup Je{fj^) < {ax,K + S)P{E, ( 52 ) 

£—^0 

Since 5 is arbitrary in (I 52 p . this concludes the proof of the upper bound. □ 

Remark 4 . 4 . With a minor adaptation of this proof (see | 35 l [ 39 ]), one could also deal with 
Dirichlet boundary conditions (that is impose jyf = 0 on d^l). Letting 

7 i = + \{Vi+r]i- if + {K - l)rilril : ryi(O) = ry2(0) = 0 , lim??* = l| , 

we would obtain as F—limit (at least for dVt of class 


Dir 


a\^KP{E, £l) + n dE) + n dE^) \i r]i = xe = I - m 

< and |L1| = ai 

^ +00 otherwise. 


Notice that in this case, by definition, 7* < crx x + Xj- Hence, on the macroscopic level, there is 
always a contact angle. Complete wetting is still possible in the form of an infinitely thin layer 
of one of the phases (this would show up in the 7*). The behavior of 7* has been recently studied 
in the physics literature | 43 j . 

Remark 4 . 5 . When K —)• +00 as e —)• 0 a refinement of the argument giving (| 51 l) gives also 
the optimal prefactor (see [IT] Prop. 6.2]). 


4.4 Study of the surface tension 

We study the asymptotics of the surface tension ax,K in the limit regimes iL —)• 1 and K —)• +00. 
We first consider the case iL —>• 1 (weak segregation). 

Proposition 4.6. There holds: 


ax,K 21 - A3 

lim , =--77. 

K^i 3 1 - A 2 

Proof. Letting x = {K — in the definition of crx,K, we have 

= int i,;i^+ xw+thl + ^ahv) 1’'!+’'L P ■ -S’'! =IfS’'! = 4 ■ 

Hence, using similar L—convergence arguments as above (see also [ 25 ]), we get that, 

I [ \Vi\^ + + V1V2 ■ lim 71 = 0 , limr/i = l,ry^ = (1 - 7 i)| • 

[Ju -°o +°o J 

Letting rji = cos 4 > for lim_oo 4 > = ^ and lim+oo </> = 0 , and using that + 6^ > 2 ab, we find 

_ roo 

lim — > 2 / |(/)'|(1 — (cos(())^(l — A^))^/^|cos(^sin(^| 

x/K - 1 J -00 


= 2 


^^^-^^(cos2(x)A 2 - cos2(x) + 1)3/2 


2 1 - A3 
31 -A 2 
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In order to reach equality it is enough to consider the unique solution of 

(/>'(! — (cos 0)^(1 — = — cos 0 sin 0, 

with (/)(0) = The solution is indeed unique and defined on R since is 

Lipschitz continuous. Moreover, it is decreasing and has the right values at ±oo. □ 

Remark 4.7. Notice that the value of limx^i exactly coincides with the one found in 

m (see also [12]). Using dal), it is moreover not hard to prove that there exists C > 0 (not 
depending on A), such that 

< ax,K < 


when K ^ 1. 

Remark 4.8. Proposition 14.61 also holds with the same proof for A = 1. The limit is then equal 
to 1. 


Eventually, we consider the strong segregation asymptotics, K —>■ +oo (notice that we recover 
the same scaling as the one predicted in the physics literature |42j 1. 


Proposition 4.9. There exist 0 < c < C such that for K > 1 and A > 0, there holds 


'^A,oo ^ 


A1/2 


+ 


{K - 1)1/4 ^1/2 


< < 0-A,oo—C 


1 


+ 


A1/2 


KI/2 7 ^ 1/4 ) ’ 


(53) 


with 


O-A.oo = inf + A^|7/2p + ^ (t/i+ r /2 - l)^ : limr?i = 0, lim??! = 1 , t/it /2 = o| 

2y/2 


= (1 + A)- 


(54) 


Proof. Let us establish the identity (fSH) . Since any admissible pair {r]i,rj 2 ) with finite energy is 
continuous and satisfies r]i{x) —)■ 1 as x —>■ +oo and r] 2 {x) ^ 1 as x ^ —oo, there exists x G R 
such that {r]i,r] 2 ){x) = 0. Using translation and symmetry, we see that 


o-A,oo = 71 + 7 a, 


with 

7a = iiif + \W ■ 7 G ([0, Too)), 7/(0) = 0, lim?/ = l| . 

The classical Modica-Mortola procedure applies to this minimization problem. For A > 0 we 
have (see e.g. Chapter 6 in M). 


7A ^ 2X^l\Y^ds ^ ^A. 


2A2 


Moreover, the minimizer is given by the formula 


T]{x) = tanh 


X 


V2X, 


for X > 0. 


(55) 
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Let us now establish the upper bound in (1531) (right inequality). If A > iL (that is 


inspired by (l55|) . we consider the prohle 
0 if X < 0, 


r]i{x) = < 


(X \ %(x) = < 

j if X > 0, 


d-x\ 

tanh —if X < (5, 

\V2X 


0 if X > (5, 

The parameter d > 0 tunes the width of the overlap between the two species (see Figure [S]). 



Figure 5: Prohle of the competitor for the upper bound 
The energy f{5) of this competitor is an upper bound for a\^K- We compute 

f{5) = Ex^xiv) = j kiP ++ 2 (^1 ^2 “ l) + (-^ “ l)hi^2 

f ^ r +00 1 ^ 

= Ih'iP + Av' 2 ? + / ^ 2^^ ■ + / 2 + i ■ 2^ 

= ^ (^^ 1^2 “ ^) = • 

Using the expressions of r/i and 7/2 and the concavity of t € R+ 1 —> tanh(t), we get 

< f{S) < ^a,oc-2 + ^ I 

Optimizing in 6, we obtain with 6 = (12A^/iL)^/^, 


-s) ds - ax,oo-^ + j^. 


(y\,K < <^\,oo 


4 /12A2\^/'^ 4(12)^/^ VA 


K 


This yields the right inequality of (I53p with c = 


^A,oo 

4(12)^/'^ 


iFV4- 


If on the contrary A <C iF we keep 771 (x) = [tanh(x/\/2)] as above and let 

'1 X < 0 

mix) = < 

•w 

The energy can then be estimated by 


1 X < 0 

l-K^/^x xe[0,iF-V2] 

0 X > iF-^/2. 


E\,Kim,V2) < 71 + A^iF^/^ + [ 

Jo 


i^-l/2 r 
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Notice now that <C K and that since tanh(x/\/2) < a;/\/2, 

r ''' [(1 - - 1] < r ''' [(1 - - 1] ’ 

JO JO 


which implies 

.K-l/2 r 


/ 


^(1 - f?2)^ + ^ 


< 


L 


X-l/2 


—x^-Kx^ = -J-K-^/^ 

2 10 


Putting all these estimates together gives the upper bound 

o'\,K < o-A,oo — . 

We now turn to the proof of the lower bound (left inequality of (I53p l. 

Let us consider an admissible pair r] = (r/i,T/ 2 ) € finite energy. Without 

loss of generality, we may assume 0 < r]i{x),r} 2 {x) < 1 for x G R. By continuity of r] and from 
the conditions at infinity, there exists xq G R such that ?7i(xo) = A? 72 (xo) = rn and a maximal 
interval Iq = [xq— 5-, xo+(^+] such that m/2 < ryi, A ?72 < 2m in Iq. We define /_ = (—oo, xq—5-) 
and /+ = (xo + 5+, +oo) and split the integration in three parts: + £'o + i?+, 

with 

Ei = f + A^|r? 2 p + - (r/i + ??2 “ l) +(-^“l)^i^ 2 i for i G { — , 0,+}. 

J 


Using the Modica-Mortola trick, we have 

E+ > 


r-\-oo rl 

/ \r}'i\^ + WKijli) ^ 2 / \/wk{s) ds. 

Jxq+ 5+ Jr]\(xQ+5+) 

Using r/i(xo + 5+) < 2m and 2 yJwK > \/2(l — s^) — V^Xio if-i/2]j we obtain, 

E+ > V 2 [ (1 - ds - > ^-2C2m-C2iL"^/^. (56) 

J2m 3 

2 C 2 A 


Similarly, 


E_ > 


-2V2m-V2K-^/‘^. 


(57) 


Now we consider the middle part. Since min ^ fj^ \ri[\, A|r; 2 l^ — ^/2, we have by Cauchy- 

Schwarz inequality, 


K? + Ari' 2 ? > 


m 


Uo 4(d+-d_)- 

On the other hand, since r/i, At /2 > m/2 in Iq, we have 

f 1 


'lo 


W + 'Ui)U(Tf-i)*| > 


Optimizing with respect to (<5+ — d-), we obtain. 


En = 


[ lVil^ + A^|r/2|^+ /" I {Vi + V2 -+ (E - l)ViV2 > 

J/o JIo ^ 


y/K — Im^ 
4A ' 


(58) 
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Gathering together (I56]) . (l^ . (f58]l . we get 

ExMri) > ^(1 + A) - 4V2m + 

Minimizing with respect to m, the minimum is reached for m = [16\/2A/(3\/-fi" ~ which 

yields, 

2 V 2 32 23/^1 A1/2 2 V 2 

> ^( 1 + )- 33/2 (ir- 1 ) 1/4 “ i^i/ 2 ' 

This establishes the left inequality of (|53|) and ends the proof of Proposition 14.91 □ 
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